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A double-acting thermoacoustic heat engine was invented recently by our group. It generally consists of 
several sequentially connected identical units. In order to apply such a system to dish-solar power gen¬ 
eration and for better absorbing concentrated solar heat, a special-shaped heater composed of a bundle of 
tubes was designed. To understand the internal operating mechanism and thermodynamic performance 
of such a solar-driven thermoacoustic generator, a test rig having only one basic unit was designed, built 
and tested. Firstly, thermoacoustic conversion performance of this test rig was theoretically predicted 
based on linear thermoacoustic model. By adjusting the external resistance and moving mass of the 
expansion motor, the acoustic filed in the THE can be regulated to fulfil better energy conversion. In 
the simulation, a highest thermal efficiency of 51.37% and a net acoustic power of 1.6 kW under a heating 
temperature of 650 °C can be obtained. Then, the experiments were conducted and compared with the 
simulation results. A discrepancy between the theoretical and experimental results was observed and 
discussed, and guidance was given to improve the experimental system. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

Solar energy is the most abundant clean energy on earth. A 
great portion of solar energy impinging on earth is thermal radia¬ 
tion. However, the efficient and reliable way to use this part of so¬ 
lar energy is very limited. Thermoacoustic heat engines (THE) can 
effectively convert solar energy to acoustic power, i.e. mechanical 
energy, based on the thermoacoustic effect, such as the Rijke and 
Sondhauss effect. THEs also have some other great advantages. 
First, a THE generally consists of several heat exchangers, regener¬ 
ator or stack and tubes, which means there are no moving compo¬ 
nents in high temperature parts, so THEs are more reliable and 
with a long life-time. Second, THEs are generally charged with inert 
gases such as helium and nitrogen, hence they are environmentally 
friendly. Third, in addition to solar energy, THEs can also use other 
low-grade energy such as waste heat. A comprehensive review of 
THEs was presented by Swift in 1988 [1]. The mechanism of THE 
can be understood based on the interaction between solid and 
gas particles. When heat is supplied through the hot end heat 
exchanger and the ambient heat exchanger is cooled, the temper¬ 
ature gradient along the regenerator or stack is established. The gas 
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particle interacts with solid materials (stack or regenerator) and 
when the temperature gradient is large enough the oscillation 
occurs, realizing the enlarging or generation of acoustic power, 
i.e. mechanical energy. With the increasing concern of high effi¬ 
cient energy conversion and the utilization of low-grade energy, 
the researches of thermoacoustic heat engines (THE) are booming 
in recent decades. 

There are two typical kinds of THEs — the standing wave ones 
and the travelling wave ones. The latter operates based on the Stir¬ 
ling Cycle which has the potential of high efficiency [2]. The first 
concept of travelling wave THE was proposed by Ceperley in 
1979 [3], and later in 1998 Yazaki et al. built a modified THE based 
on Ceperley’s work 4]. But their THEs had a serious problem of 
large viscous loss in the regenerator or stack. A breakthrough of 
travelling wave THE was made in 1999 when Backhaus and Swift 
[5] developed a travelling wave THE with a thermal efficiency of 
30%. In 2005, Luo et al. [6] proposed an energy-focused travelling 
wave THE which contained a tapered resonance tube to reduce 
the nonlinear dissipation in the resonance tube and the pressure 
ratio reached 1.4. In 2010, Kees de Blok proposed a novel 4-stage 
travelling wave THE 7]. The multi-stage configuration is beneficial 
for even higher energy density and larger power capacity. There 
are also some real applications for THEs. Since 1994, a project 
focusing on the application to natural gas liquefier driven by THE 
has been conducted [8 . A thermoacoustic generator for 
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Nomenclature 



Abbreviation 

Greek symbols 

THE 

thermoacoustic heat engine 

y 

specific heat ratio 

1-D 

one dimension 

3 

penetration depth, m 



6 

phase delay, ° 

Symbols 


( 

correction factor of channel solid thermal parameters 

a 

speed of sound, m/s 

P 

gas density, kg/m 3 

A 

area, m 2 

G 

Prandtl number 

B 

magnetic field, T 

T 

transduction coefficient 

c 

heat capacity per unit mass, J/kg K 

CO 

angular frequency, s _1 

D 

diameter, m 

Re[] 

real part of 

E 

acoustic power, W 

Im[] 

imaginary part of 

f 

frequency, Hz 

r >sj 

conjugation of a complex quantity 

i 

spatially averaged diffusion function 

1 1 

amplitude of a complex quantity 

H 

total energy power, W 



i 

V-l, imaginary unit 

Subscript 


I 

input current, A 

1 

first order 

k 

thermal conductivity, W/m K 

2 

second order 

I< 

spring constant, N/m 

e 

electrical 

l 

length, m 

m 

mean 

L 

electrical inductance, H 

m 

mechanical 

M 

mass of piston, kg 

P 

isobaric 

P 

pressure amplitude, MPa 

s 

solid 

P 

working pressure, MPa 

K 

thermal 

Q 

heat input per unit length, W/m 

V 

viscous 

Q 

heat transferred, W 

com 

compressor 

R 

electrical resistance, C2 

gen 

expansion motor 

R 

mechanical resistance, N s/m 

in 

input 

T 

temperature, I< 

in 

inlet 

U 

volume flow rate, m 3 /s 

net 

net 

V 

input voltage, V 

out 

outlet 

z 

impedance, Pa s/m 3 




applications such as LED lightening or charging battery driven by 
waste heat from stove has been constructed recently [9]. However, 
the above-mentioned THEs have a common problem of having a 
long resonant tube, which has limited the further application of 
THEs. Thus, in recent time our research group has presented a con¬ 
cept of double-acting multi-cylinder travelling wave THE [10] to 
solve the problem. It is schematically shown in Fig. 1. This engine 
mainly consists of several sequentially connected basic units in a 
closed loop. Each basic unit is composed of main components of 
a THE sandwiched by the compression and expansion pistons. 
“Double-acting” means one piston works as compression piston 
to the previous THE unit and simultaneously as expansion piston 
to the last unit. This novel configuration not only has the tradi¬ 
tional advantages of THEs but also some other advantages as stated 
below. This double-acting THE is capable of producing larger 
power because of its multi-cylinder configuration. In addition, sev¬ 
eral symmetric units cascading in one loop is equivalent to several 
THEs working together, hence the configuration is relatively com¬ 
pact. Taking advantages of high inertia of the pistons this configu¬ 
ration also eliminates the long resonant tube. Recently, we are 
considering constructing a double-acting THE for dish-solar power 
generation, aiming at a 5 kW acoustic power output which is for 
used to drive electric expansion motors. In order to effectively ab¬ 
sorb the concentrated solar heat, a special-shaped heater com¬ 
posed of a bundle of tubes was designed. Different from the 
generally-used shell-and-tube and plate-fin heat exchangers, this 
heater has much longer and less tubes and the shape is irregular 
which will result in problems of jet-flow and turbulence. The 
whole system is complicated, but it consists of several symmetric 


units which mean each basic unit should have the same perfor¬ 
mance. Hence in order to study the advantages as well as draw¬ 
backs, and find ways to improve it, a test rig having only one 
basic unit which is representative to the system was designed, 
built and tested. Firstly, a numerical analysis was carried out to 
study the thermal conversion performance. Then experiments 


6 



1- main ambient heat exchanger; 2-regenerator; 3-heater; 

4-thermal buffer tube; 5-secondary ambient heat exchanger 
6-connection tubes; 7- motor 


Fig. 1. Schematic of double-acting THE. 
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were conducted and compared with the simulation results. Finally, 
conclusions were made. 

2. Test rig 

Fig. 2 is an assembly drawing of a 6 x 1 kW double-acting THE 
used for dish solar power generation. This 6x1 kW double-acting 
THE consists of two 3x1 kW double-acting THEs, and the two 
sub-systems are carefully arranged to minimize the vibration. 
The top part of the system is a heater composed of a bundle of 
tubes. This configuration is beneficial for absorbing concentrated 
solar heat. 

Since the configuration is complicated, a test rig having only 
one basic unit is built to understand the internal operating mech¬ 
anism. The schematic diagram of the test unit is shown in Fig. 3. 
The test unit consists of a compression motor, a THE and an expan¬ 
sion motor. The THE consists of a main ambient heat exchanger, 
regenerator, heater, thermal buffer tube and secondary ambient 
heat exchanger. The compression motor and expansion motor are 
both of dual-opposed configurations. The two pistons of each 
motor move in the opposite direction and since their same config¬ 
urations and same mass the vibration is reduced based on the 
conservation of momentum. The compression motor and expan¬ 
sion motor are working simultaneously to modulate the acoustic 
field. The acoustic power is input from the compression motor, 
amplified in the tested THE unit and then consumed in the expan¬ 
sion motor. 

Consequently the net acoustic power E net generated in the 
tested THE unit is: 

Enet = Egen ~ E com (1) 

where E gen is the acoustic power consumed in expansion motor and 
E C om is the acoustic power input from the compression motor. 

And the thermal efficiency of the tested THE unit is: 

f? = ^xl00% (2) 

Ccin 


where Qin is the heat input in the heater. 

The impedance of motor can be described as: 



Fig. 2. A 6 x 1 kW double-acting THE used for dish solar power. 


1 7 



1- compression motor; 2- main ambient heat exchanger; 

3- regenerator; 4- heater; 5- thermal buffer tube; 

6- secondary ambient heat exchanger; 7- expansion motor 

Fig. 3. Schematic of a basic unit of a double-acting THE used for dish solar power 
generation. 

where R e is the electrical resistance, including the internal resis¬ 
tance, i.e. the resistance of the coil, and the external resistance. 

The distribution of acoustic filed, which is determine to the 
acoustic power generated, can be mathematically reflected by 
the distribution of pressure amplitude P\ and volume flow rate 
U i, i.e. the product of oscillation velocity and gas flow area. Accord¬ 
ing to Eq. (3), by adjusting the impedance of the expansion motor, 
we can obtain different acoustic distributions, hence the perfor¬ 
mance of the system under different acoustic conditions could be 
studied. Although many parameters of the motor can be adjusted 
to realize the adjutancy of impedance, the change of the external 
resistance R and mass of the piston M is most convenient. Conse¬ 
quently, in the following studies, we adjust these two parameters 
to establish different acoustic conditions to study the acoustic 
power generation in the tested THE unit. 

Table 1 gives the dimensions of the tested THE unit. Table 2 
gives the parameters of the compression motor and expansion mo¬ 
tor. These dimensions are coincident in the 6x1 kW system, the 
numerical model and the test rig. Since the compression motor 
and the expansion motor each have two same pistons, only one 
single piston’s parameters are given. 

3. Simulation analysis 

3.1. Simulation model 

Based on the linear thermoacoustic theory, the simulation is 
conducted using the DeltaEC 6.2 program [11]. This program pro¬ 
vides many segments including pipe, regenerator, heat exchanger, 
etc. for users to simulate different components of a thermoacoustic 
system. The program integrates the governing equations in one 
dimension (1-D). The geometry parameters are given by the user. 
And the thermophysical properties of gas and solid are given in 
the DeltaEC’s library. A shooting method is introduced in the pro¬ 
gram to satisfy various boundary conditions set by the users. 

3.1.1. Geometric model 

The simulation model is based on the construction of the basic 
unit, but some simplification has been made, as seen in Fig. 4. 
Firstly, the real compression motor is dual-opposed, and the 
electric circuits of the two pistons are in parallel. In the simulation, 
we adopt only one piston which is equivalent to the two paralleled 
pistons. Hence, the cross-sectional area of this equivalent piston is 
the sum of the two pistons’ areas. And other parameters of the 
equivalent piston equal to the sum of the corresponding 
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Table 1 

Dimensions of tested THE unit. 



D (mm) 

l (mm) 

Others 

Main ambient heat exchanger 

55 

52 

Shell-and-tube heat exchanger 241 tubes with 1.5 mm inner diameter 

Regenerator 

55 

40 

150 mesh stainless steel screen 

Heater 

2 

210 

19 tubes 

Thermal buffer tube 

38 

120 

/ 

Secondary ambient heat exchanger 

38 

30 

Shell-and-tube heat exchanger 121 tubes with 1.5 mm inner diameter 


Table 2 

Parameters of compressor and generator (single piston). 



D (mm) 

M (kg) 

K (1<N m” 1 ) 

R m (Nsrn t ) 

R e (Q) 

Bl (NA' 1 ) 

Compressor 

60 

2.08 

180 

16 

3.72 

97.7 

Generator 

50 

1.1 

167.5 

20 

3.5 

95 


Working Conditions 


y=80 Hz 
Pm=4 MPa 



1- equivalent compression motor; 2- main ambient heat exchanger; 

3- regenerator; 4- heater; 5- thermal buffer tube; 

6- secondary ambient heat exchanger; 7- equivalent expansion motor 

Fig. 4. Geometric model. 


parameters of the two paralleled pistons. For the expansion motor, 
we adopt the same treatment. Secondly, we simplify the special¬ 
shaped heater into a shell-and-tube heat exchanger with an inner 
diameter of 55 mm and with 19 tubes whose inner diameter is 
2 mm. 

3A 2. Governing equations 

For the components charged with gas, based on the linear ther¬ 
moacoustic theory [12], the governing equations are: 

dlh iapnlh 

dx 1 -/„ A w 


where oo is angular frequency and co = 27 if, here the resonant fre¬ 
quency/is set as 80 Hz./ V and / K are viscous function and thermal 
function. Other symbols are defined in Nomenclature in detail. 

At the inlet of one segment, p a , Lfi, T m equal to the outlet param¬ 
eters of the previous segment. P m is set to be 4.0 MPa. A, A s are gi¬ 
ven uniquely in different segments which can be calculated from 
the dimensions given in 'ables 1 and 2. q equals to 0 in the most 
of the segments except for the heat exchangers. In the heat 
exchangers, q are solved to meet the targets of the specific solid 
temperatures in the heat exchangers. 

With the assumed values of p A and Lfi, the acoustic power E can 
be calculated: 


dU A 

~dx 


x 


icoA 

W m L 

Ui dT 


1 + 


m 


T m dx 


(y-i)/ K 
l + £ 


Pi + 


(f K ~fv) 


(1 -/„)(! -<x)(l+{) 


(5) 


dH- 
~dx 

where 


= Q 


( 6 ) 


Hi =1 Re 


+ 


Pil/i 1 


Pm C pl^l 


f K ~fi 




dT 


m 


2Aco(\ - g)\\ - f v 


dx 


Im 


fv + 


(fK-fo)V+tfo/f K ) 

(!+{)(!-(7) 


— (Ak -(- A s k s ) 


dT 


m 


dx 


E = ^Re\pA 


1 

2 


RejpiLTi] 


1 

2 


|p 1 1 |f/i | cos 6 



where 9 is the phase difference between pi and Lfi. 

For the pistons, the VESPEAKER model is used. The governing 
equations are: 


II — II C ° ' 1 ^ K A n 

1/1 Uhm p m a 2 1 +£ 2 APl m 

(8) 

Pi, out — Pi,in T ^1 — Z m U\ 

(9) 

Vl = Ze/i - Tl/! 

(10) 

Z e = R e + icol 

(11) 

= R m />i 2 + i(a>M - K/co)/A 2 

(12) 
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t = -t ' = Bl/A (13) 

where Lf 1>in is the volume flow rate at the inlet of the VESPEAKER, 
Pi.out. Pi,in are the pressure amplitudes at the outlet and inlet of 
the VESPEAKER, respectively. We balance L with an external electri¬ 
cal compliance, hence L is 0 here. For compression motor, V\ is 
guessed to meet the target of the displacement amplitude of the 
expansion motor piston. For expansion motor, V\ is 0. 

3.1.3. Boundary conditions 

The beginning of the simulation model is the back cavity of the 
compression motor; hence the inlet Eh is 0. In the main and sec¬ 
ondary ambient heat exchangers, the solid temperature is set as 
20 °C. The solid temperature of the heater is 650 °C. We consider 
that the piston of the expansion motor is a thermally insulated 
model, so H 2 = E at the inlet of the expansion motor. The piston dis¬ 
placement amplitude of the expansion motor is set as 6 mm to pre¬ 
vent the piston from hitting the cylinder due to the limitation of 
back cavity volume. The end of the simulation model is the back 
cavity of the expansion motor, the outlet Eli is 0. 

3.1.4. Shooting method 

In the DeltaEC program, a shooting method is introduced in the 
program to satisfy various boundary conditions set by the users. 

At the inlet of the simulation model, we set the values of 
pressure amplitude p gas temperature T m , the input voltage V\ 
of the compression motor, and the phase angle of the input voltage 
of the compression motor as “guessed parameters” and give them 
initial values. By integrating the governing equations, the outlet 
parameters of the simulation model can be obtained. If the outlet 
parameters do not meet the targets, the guessed parameters will 
be modified and given to the inlet of the model again until the out¬ 
let parameters meet the target parameters. 

3.2. Simulation results and discussion 

3.2.1. The influence of external resistance and moving mass 

The distribution of acoustic field is not only influenced by the 
geometry parameters of the THE but also largely affected by the 
impedance of the expansion motor. Based on the above analysis, 
we can regulate the acoustic filed by adjusting the external resis¬ 
tance and the moving mass of the expansion motor. In this section, 
we mainly study the performance of the basic unit when the exter¬ 
nal resistance and the moving mass are changed. Note that the 
following-mentioned “external resistance” is the resistance which 
has already been transformed to the resistance needed in the 
paralleled electric circuit in the experiments. The following-men¬ 
tioned moving mass is the mass of a single piston. In the following 
simulations, the external resistance varies from 5.25 Q to 110 Q. 
When the external resistance is smaller than 5.25 Q, a positive 
thermal efficiency and net acoustic power are unable to be ob¬ 
tained. Three specific moving masses are studied, i.e., 1.0 kg, 
1.1 kg and 1.2 kg. 

The phase difference of volume flow rates at the interfaces of 
the compressor piston and the expansion motor piston is dominant 
because it can reflects the acoustic filed distribution in the tested 
THE unit. The distribution of acoustic field will further influence 
the performance of the basic unit, especially the energy conversion. 
Fig. 5 illustrates the curves of the phase difference of volume flow 
rates at the interfaces of the compression motor piston and the 
expansion motor piston. As seen in the figure, both the moving 
mass and the external resistance have a significant influence on 
the phase difference of the volume flow rates. When the moving 
masses are 1.0 kg, 1.1 kg and 1.2 kg, the phase differences range 
from -100.74° to -139.09°, from -104.88° to -152.72° and from 
-108.93° to -160.73°, respectively. 



External Resistance (&) 


Fig. 5. External resistance versus phase delay of volume flow rate between pistons. 

The net acoustic power generated and the thermal efficiency of 
the tested THE unit are the most important two parameters that 
imply the performance of the basic unit. An extremum of net 
acoustic power and thermal efficiency is obtained with suitable 
external resistance, which means the best acoustic field distribu¬ 
tion has been obtained. When the external resistance is smaller 
than 10 Q, the acoustic power generated in the tested THE unit 
abruptly jumps to the maximum value with only a bit increase of 
the external resistance, as shown in Fig. 6. With the further in¬ 
crease of the external resistance, the net acoustic power decreases 
smoothly. The external resistance versus thermal efficiency dia¬ 
gram has similar trend to that of the net acoustic power, as shown 
in Fig. 7. When the external resistance is smaller, the thermal effi¬ 
ciency increases dramatically to a maximum value with the in¬ 
crease of external resistance. Then the thermal efficiency has a 
slight decrease. However, the best acoustic field distribution is dif¬ 
ferent to acoustic power and thermal efficiency. Hence the choice 
of external resistance should be compromised. 

Similarly, there should also be a most suitable moving mass 
which will result in the extremum of net acoustic power and ther¬ 
mal efficiency. However, of the three moving mass we choose, i.e. 
1.0 kg, 1.1 kg and 1.2 kg, the most suitable moving mass has not 
obtained yet. 

Among the current moving mass and external resistance, the 
maximum acoustic powers generated are 1566.6 W, 1592.5 W 
and 1600.5 W with an external resistance of 10 Q when the mov¬ 
ing masses are 1.0 kg, 1.1 kg and 1.2 kg respectively. When the 



External Resistance ( Q ) 

Fig. 6. External resistance versus acoustic power with different moving mass. 
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External Resistance ( Q) 


Fig. 7. External resistance versus thermal efficiency with different moving mass. 

moving masses are 1.0 kg, 1.1 kg and 1.2 kg, the highest thermal 
efficiency of 47.72%, 49.81% and 51.37% are achieved with 29 Q, 
32 Q and 42 Q external resistances respectively. 

In general, for a higher net acoustic power and thermal effi¬ 
ciency, a relatively smaller external resistance and a relatively 
higher moving mass is preferred. 



Diameter of heater tubes (mm) 

Fig. 9. Diameter of heater tubes versus acoustic power and thermal efficiency. 

increase of heat transfer area will contribute to the increase of gen¬ 
erated net acoustic power and thermal efficiency. However, as the 
diameter is larger than 3 mm, the surface viscous loss outweighs 
the heat transfer area, which results in the slight decrease of the 
net acoustic power and thermal efficiency. When the inner diame¬ 
ter of the tubes is 3.1 mm, the highest thermal efficiency of 51.19% 
and the maximum net acoustic power generated of 846.28 W have 
been achieved. 


3.2.2. The influence of heater 

In this section, we mainly discuss the influences of the heater. 
The moving mass of the expansion motor is 1.1 kg and the external 
resistance is 32 Q. 

As mentioned above, the heater consists of 19 tubes with inner 
diameter of 2 mm. The influences of the number of tubes and tube 
diameters are very similar as they both determine the heat transfer 
area between gas and solid. Therefore, we only study the influence 
of the tube’s inner diameter on the performance of the basic unit in 
this section. The inner diameter of the tube varies from 1 mm to 
6 mm while the number of the tubes is fixed as 19. In the threshold 
of diameter, the phase difference varies from -147.39° to 
-128.76°, as shown in Fig. 8. The diameter of the heater tubes 
has a strong influence on the phase difference when the diameter 
is relatively small, i.e., smaller than 3 mm. 

Fig. 9 demonstrates the influence of the diameter on net acous¬ 
tic power generated in the tested THE unit and the thermal effi¬ 
ciency. When increasing the tube diameter, the heat transfer area 
increases. Hence the enhancement of heat transfer caused by the 



4. Experiment analysis 

In this section, preliminary experiments were conducted and 
discussed. The experimental results are compared with the calcu¬ 
lated results. Although there are still some issues left and the 
experimental results have much room for improvement, the funda¬ 
mental trend of the experimental results is in reasonable agree¬ 
ment with the calculated results. 

4.1. Experimental system 

The implementation of the experimental system is shown in 
Fig. 10. Six silicon carbide heating rods are used for heat input. 
The rated power of each single silicon carbide heating rod is 
800 W. Polycrystalline mullite is used for thermal insulation. The 
system is charged with 4.0 MPa helium and the resonant frequency 
is 80 Hz, which is the same as the simulation. A sheathed thermo¬ 
couple is located at the surface of the heater on the regenerator 
side to measure the heater solid temperature. The heater solid 
temperature is set as 650 °C approximately. The temperature of 



Fig. 8. Diameter of heater tubes versus phase difference of the volume flow rate. 


Fig. 10. Experimental system with heat insulator. 
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the gas in the heating chamber surrounded by the Polycrystalline 
mullite and the solid temperature of the hot end of the regenerator 
are measured to estimate the heat leakage. The pressures of the 
back and front cavity of the compressor and the expansion motor 
are measured. 

4.2. Experimental results and discussion 

A preliminary experiment is conducted and the results are com¬ 
pared with the simulation results in this section. Limited by the 
current of the compression motor, the external resistance can only 
vary from 40 Q to 110 Q. 

Figs. 11 and 12 give the dependence of net acoustic power gen¬ 
erated and thermal efficiency on the external resistance, respec¬ 
tively. The experiment results and the simulation results have a 
similar trend. However, there is a large deviation between the 
experiment results and simulation results. Especially for the ther¬ 
mal efficiency, in the experiment the thermal efficiency is quite 
small and less than 10%. It indicates that the heat input through 
the heater might have an even larger deviation between experi¬ 
ments and simulations. 

Fig. 13 shows the volume flow rate amplitudes at the compres¬ 
sor outlet and at the expansion motor inlet. Fig. 14 shows the pres¬ 
sure amplitudes at these two positions. These two positions can 
also be approximately considered as the inlet and outlet of the 
tested THE unit. Since the piston displacement of the expansion 
motor is set as 6 mm, the volume flow rate at the expansion motor 
side is independent on external resistance. In Fig. 13 we could see 
that the experimental results of the volume flow rate amplitude is 
close to the simulation results. However, in Fig. 14, it is noticeable 
that the pressure amplitude drop from the compressor side to 
expansion motor side is much larger in the experiments than in 
the simulations. That’s mainly because in the simulations, we sim¬ 
plify the special-shaped heater as a shell-and-tube heat exchanger 
so that the losses are much smaller. 

There are some possible reasons that caused the large deviation 
between the experimental results and simulation results. Some 
possible ways to reduce the deviation are discussed here. 

1. The structure of the heater has caused a very large flow loss and 
pressure drop. According to Fig. 14, the pressure drop in the 
experiment is much larger than that in the simulation. That 
means the pressure drop in the heater is much larger than 
expected. The sudden change of cross-sectional areas from a 
large tube to several small tubes will cause significant turbu¬ 
lence in the heater. Hence the acoustic power is terribly con¬ 
sumed in the heater. However, in the simulation, this loss is 






ignored to a large extent. By adding some rectifiers at the outlet 
of the heater may find help to suppress the turbulence to some 
extent. 

2. The large deviation is mainly because of the heat loss and low 
heat absorption efficiency of the heating system. The total input 
electric power of the six silicon carbide heating rods is consid¬ 
ered as the amount of the heat input into the system. However, 
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the practical heat used by the test rig is much less though it is 
difficult to determine how much energy is actually used by the 
system. First, about 30% of the length of the rods is not insulated 
in the Polycrystalline mullite for circuit connection. That means 
at least 30% of the heat radiation from the rods is wasted. Sec¬ 
ond, the radiation efficiency is limited by the confinement of 
the surface temperature of the rods, the radiation angle and 
the surface gray level of the heater. This results in a large tem¬ 
perature difference between gas and solid, thereby resulting in 
a large amount of heat wasted, more than 50%. Fig. 15 shows an 
infrared photo of the basic unit when it is heated. Fig. 16 shows 
the temperature variation with the external resistance. The 
temperature of the heater on the regenerator side is set as 
650 °C approximately. The temperature of the gas in the heating 
chamber and the solid temperature of the hot end of the regen¬ 
erator are also given in the figure. From Fig. 16 we can learn that 
the heat leak is very serious. The gas temperature in the heating 
chamber varies from 810.6 to 889.3 °C while the heater solid 
temperature is only 650 °C. Only 1949 W heat input is needed 
to maintain the solid temperature of the heater at 650 °C when 
the test rig is not working. While when the test rig is working, 
the heat input is 4341 W with an external resistance of 40 CL 
Based on the above analysis, about 50% of the heat is lost 
because of the heat leak. Third, the sudden change of cross-sec¬ 
tional areas from the heater to the thermal buffer tube will 
cause a jet-flow in the thermal buffer tube. This jet flow will 
cause heat losses through the thermal buffer tube. All these fac¬ 
tors causes the serious heat leakage. Since the heat leak is so 
serious, much more heat should be input to maintain the tem¬ 
perature of heater. The thermal efficiency will have a terribly 
decrease consequently. 



External Resistance ( Q) 

Fig. 16. External resistance versus temperature. 


having a 100-150 °C temperature difference. With a much 
smaller temperature gradient in the regenerator, the acoustic 
power generated is smaller. 

Generally speaking, the simulation is based on the ideal condi¬ 
tions to a large extent. Although there is a large deviation between 
the calculations and the experiments, the trends are coincident. 
Moreover, there is still large room for the experimental results 
improvement. Although some issues are left, the system shows 
its potential for application in solar power generation. 


To solve this problem, some methods can be adopted. First, a 
better heating method such as heating through conduction should 
be considered. Second, a better insulation is necessary. Third, some 
meshes can be located at the inlet of the thermal buffer tube to 
suppress the jet-flow. Also the diameter of the heater tubes should 
be optimized. 

3. Fig. 16 also reveals that there is a large temperature difference 
between the heater solid and regenerator solid. We set the solid 
temperature of the heater as 650 °C. In the simulations, the gas 
temperature at the high temperature side of the regenerator is 
near 650 °C since the heat transfer is relatively ideal. Flowever, 
in the experiments, due to the abrupt change of the cross- 
sectional area from the regenerator to heater, the flow at this 
interface is quite complicated. The gas temperature at the high 
temperature side of the regenerator is much lower than 650 °C, 



Fig. 15. The infrared photo of the basic unit when heating. 


5. Conclusions 

In this paper, a 6 x 1 kW double-acting TF1E is considered for a 
further application to drive generators to fulfil the dish-solar 
power generation. A basic unit which is representative to the 
whole system is numerically and experimentally studied. Some 
conclusions are drawn as follows. 

1. The performance of the test rig was theoretically predicted. By 
adjusting the external resistance and moving mass of the 
expansion motor, the acoustic filed in the THE can be regulated 
to fulfil better energy conversion. A highest thermal efficiency 
of 51.37% and a net acoustic power of 1.6 kW under a heating 
temperature of 650 °C have been obtained with suitable exter¬ 
nal resistance and moving mass. 

2. The heater parameter has a large influence on the performance 
of the basic unit. 

3. Experiments have been conducted and the experimental results 
have been compared with the calculated results. Though some 
issues are still left, the system show its potential for an applica¬ 
tion for electric generation. Some improvements can be made to 
enhance the experimental results, such as improving the insula¬ 
tion, finding better ways to heat and adding some rectifier or 
meshes. 

The further optimization and improvement of the heater will be 
carried out in near future to improve the thermal performance of 
the test rig. 
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